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Abstract 

The Hubbard Hamiltonian and its variants/generalizations continue to domi- 
nate the theoretical modelling of important problems such as high temperature 
superconductivity. In this note we identify the set of relevant operators for the 
Hubbard Hamiltonian with a magnetic field term. 



*Permanent address: Department of Physics, University of Peshawar, Peshawar, NWFP, Pakistan. 
^Institute of Applied Physics, Tsukuba University, Tsukuba, Ibaraki 305, Japan. 
^Institute of Information Sciences and Electronic, University of Tsukuba, Tsukuba, Ibaraki 305, 
Japan. 

1 



I. INTRODUCTION 



The Hubbard Hamiltonian [HH] and its extensions dominate the study of strongly 
correlated electrons systems and the insulator metal transition JTJ. One of the attractive 
feature of the Hubbard Model is its simplicity. It is well known that in the HH the band 
electrons interact via a two-body repulsive Coulomb interaction; there are no phonons in 
this model and neither in general are attractive interactions incorporated. With these points 
in mind it is not surprising that the HH was mainly used to study magnetism. In contrast 
superconductivity was understood mainly in light of the BCS theory, namely as an instability 
of the vacuum [ground-state] arising from effectively attractive interactions between electron 
and phonons. However Anderson suggested that the superconductivity in high T c material 
could arise from purely repulsive interaction. The rationale of this suggestion is grounded 
in the observation that superconductivity in such materials arises from the doping of an 
otherwise insulating state. Thus following this suggestion the electronic properties in such 
a high T c superconductor material close to a insulator-metal transition must be considered. 
In particular the one-dimensional HH is considered to be the most simple model which can 
account for the main properties of strongly correlated electron systems including the metal- 
insulator transition. Long range anti-ferromagnetic order at half-filling has been reported in 
the numerical studies of this model HH . Away from half-filling this model has been studied 
in 1,1. 

The Maximal Entropy Principle [MEP] is a useful tool to get the dynamical and 
thermodynamical descriptions. The main advantage of this formalism is to provide a definite 
prescription to determine the complete set of operators [i.e. relevant operators] related to the 
problem under considerations. An attractive feature of the relevant operators is that they 
are group theory based and hence once a Hamiltonian for a system can be written down, 
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the task of identifying the relevant operators can proceed in principle. Relevant operators 
for a given physical system along with the Hamiltonian in essence describe the essential bare 
bones of the physical system. In a series of papers 0HH> ^he generalized time- dependent 
Jaynes-Cummings Hamiltonian in the context of Maximum Entropy Principle [MEP] and 
group theory based methods fH| was studied. In particular, in the MEP formalism 
was used to solve time-dependent N-level systems. A set of generalized Bloch equations, in 
terms of relevant operators was obtained and as an example the N = 2 case was solved. It 
was thus demonstrated in [|7| that the dynamics and thermodynamics of a two-level system 
coupled to a classical field can be fully described in the framework of MEP and group theory 
based methods. Further in || a time-dependent generalization of the JCM was studied 
and by showing that the initial conditions of the operators are determined by the MEP 
density matrix the authors were able to demonstrate that inclusion of temperature turns 
the problem into a thermodynamical one. An exact solution was also presented in the time 
independent case. Finally in || more detailed analysis of the three set of relevant operators 
was given. These set of operators are related to each other by isomorphisms which allowed 
the authors to consider the case of mixed initial conditions. The mean values of the field's 
population, correlation functions and nth-order coherence functions are of interest and useful 
in several applications. The MEP formalism allows us to describe a Hamiltonian system in 
terms of those, and only those, quantum operators relevant to the problem at hand. Thus, 
this formalism is suitable to study the Hamiltonian given in In [|||9j the population 

of each level and not their difference is considered therefore the resulting Hamiltonian is 
called a generalized time-dependent JCH. Recently the relevant operators for the generalized 



time-dependent m-photon Jaynes-Cummings Hamiltonian were determined in [II]. 

The HH for different band-fillings is studied in the context of MEP by Aliga and Proto 



12 . The HH with a magnetic field was considered by Alam and Proto [TL3| using MEP 
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techniques. In the present note we incorporate a magnetic field term in HH and identify the 
relevant operators. The set of relevant operators and their evolution equations without the 
magnetic field term considered by Aliga and Proto has also been independently checked by 
us. Moreover by neglecting the magnetic field term we easily recover the case considered 
by Aliga and Proto which provides a check on our calculations. It is interesting to note 
that Essler et al. have suggested an extended Hubbard model which contains the 

t-J model as a special case. In fact this model is a mixture of the Hubbard and the t-J 



model. The model of Essler et al. |L4[] contains a magnetic field term. On a one-dimensional 
lattice Essler et al. |14| present an exact solution to their model via Bethe ansatz. It is 
further claimed that by using ^-pairing mechanism one can construct eigenstates of the 
Hamiltonian with off-diagonal long-range order and that in the attractive case the exact 
ground state is superconducting in any numbers of dimensions. The model of Essler et 



al. |TJ,|T5[ is motivated by high-T c superconductivity and is expected to describe a system 
of strongly correlated electrons. The model of Essler et al. |14]] possesses a huge symmetry 
group [for example it has eight supersymmetries] and it would be interesting to obtain the 
set of relevant operators corresponding to it. 

The main purpose of this paper is to answer the question: Can we identify a set of 
relevant operators for the Hubbard Hamiltonian including a magnetic field term? The aims 
of this short note is to give such a set and the evolution equations for it. The layout of 
this paper is as follows. Section two contains discussion and definitions relevant to Hubbard 
Hamiltonian in the context of mean field method. In section three we recall some well-known 
results of the group theory based MEP formalism. In section four we give the relevant 
operators and the evolution equations for their expectation values in the context of the 
Hubbard Hamiltonian without and with an a magnetic field term present. Conclusions are 
given in the last section. 
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II. HUBBARD HAMILTONIAN AND THE MEAN FIELD METHOD 



The Hubbard Hamiltonian can be written as 

<i,j>, a i 

T is the hopping parameter between the nearest neighbours, c\ a creates an electron with spin 
a at site i, Ci a destroys an electron with spin a at site i, U is the on-site Coulomb interaction 
and fiirj = c\ a Ci a is the number operator for spin a at site i. 
A modified Hubbard Hamiltonian, 

H 1 = -r £ &%* + UY^fci-h[nil-k\, (2) 

<i,j>, cr i 

is also used by some authors 0,0. One may rewrite Hi in terms of H 

Hi = -r £ Lcja + UY^[^n ii -^(h^+h il ) + -} 

<i,j>, cr i 4 

= ^-?Eh+^l + ^- ( 3 ) 

i 

N in Eq. |3| is the number of sites. It is important to note that the Hamiltonians given in 
Eq. [I] and Eq. |^ cannot be considered as equivalent even when they lead to the same set of 
the relevant operators, since the g matrix [see Eq. |22| below] associated with the Hamiltonian 
in Eq. [I] is different from that which corresponds to the Hamiltonian in Eq. |3|. 

In order to get solvable model, in this note, we resort to the mean-field method. Our 
main approximation is to replace the product of operators in the hopping term by averages 
according the rule 

AB =< A > B + A < B > - < A >< B > . (4) 
It is important to note that in contrast to [TH] the mean-field approximation in our case, like 



17 1, has been applied to the hopping term. In [16] the mean- field approximation is applied 



to the Coulomb term. In our approximation the Hamiltonian can be written in site-diagonal 
form, the sites being coupled only by the mean-field parameter A a , for the definition of A CT , 
see Eq. |7| below. 

In order to apply the above rule, viz Eq. |] to the hopping term we rewrite the latter 

as 

<i,j>, a 

= \piu + ^jali^io- "I" c ja] ^ia f^jai (5) 

Applying the definition of the averaging procedure, viz Eq. ^ to the hopping term 
written as in Eq. |5| we obtain 

5Z ^icrCjcr 5Z ^ia^jc ^jcr^icri 
<i,j>, cr cr 

w K < 4 + 4 > +A CT < c ia + c jff > 

(7 

-| < A a > \ 2 I - h ia - h ja , (6) 
where we have used the definitions 

k = < 4+4 >> 

A; = < c iCT + c jff > . (7) 
Using the reduction given in Eq. |6] we may write the Hubbard Hamiltonian in Eq. [L] as 
H = T ' £ ^ + %v - t' S A a[c iCT + c i(T ] - r 53 A ^[4 + 4l 

cr cr cr 

+r'5Zl < A - > l^I + U^n^hn. (8) 

Next we want to write the HH in terms of one site only to this end we observe that 
in one dimension [1-d] each site has two nearest neighbours, in 2-d each site as 4 nearest 
neighbours and so on. Denoting the number of nearest neighbours by m we define 
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def ' 

t = mr , 

[A|2 g |A T l a +|A;| a (g) 

Using these definitions we may write the one-site equivalent of Eq. |^ 

Hi = T fl ia - T J2 A «Ci<7 ~ T J2 A *Av 
a a a 

+t\A\ 2 I + Uh i1 h il . (10) 
Eq. [T_rj allows us to rewrite the Hamiltonian in Eq. |3] in the form 

Hu = (r - — ) flier ~ r J2 A CT c ia - r K^t 

+ (r\A\ 2 + ^)i + Un^n n . (11) 

It is convenient to introduce a compact notation 

Hi = ahi - rxi + jl + Uf h (12) 

where we have defined 

~ def ~ 

rii = riii + n ib 

Xi = A*cJ T + Alcl + A T c iT + A ; ci|, 

fj = f h^hii. (13) 

hi is the number of electrons at site i, Xi is the mean field hopping interaction between 
neighbouring sites and measures the double occupancy probability or simply the number 
of pairs at the site i. 

The Hamiltonian in Eq. |Ij is a special case of the Hamiltonian form in Eq. |12| with 
the identifications a = r and 7 = t|A| 2 . If we set a = r — ^ and 7 = r|A| 2 + j in Eq. [12] 
we recover the Hamiltonian form given in Eq. ITT]. 



The magnetic field can be readily accommodated by adding the term h(n^ — h^) to 
the Hamiltonian form in Eq. [12], viz, 

Hi = ahi - TXi + jl + Ufi + hN i: (14) 

where we have defined 

Ni = ~ ( 15 ) 

fij and iVj respectively represent the symmetric and antisymmetric sums of the number 
operators for both types of spin. We define |A| 2 in analogy with |A| 2 as 



A 2 — f 


IAtI 2 -!^! 2 


2 


|A T | 2 = 


|A| 2 + A| 2 , 


|A,| 2 = 


|A| 2 - A| 2 , 


|AJ 2 in 


terms of A 2 and 



(16) 



III. OUTLINE OF THE MEP FORMALISM 



It is instructive to summarize the principal concepts of the MEP [f7|-|9| Jl8| , ^9|] ■ A sum- 



mary of MEP formalism has been given in [TT| . Here we again outline it for the benefit of 



the readers not familiar with [11 



Given the expectation values <Oj> of the operators Oj, the statistical operator p(t) 
is defined by 

/3(t)=expf-A I- S A Aj > ( 17 ) 



S 



where L is a natural number or infinity, and the L+l Lagrange multipliers A.,-, are determined 
to fulfill the set of constraints 

<6 J >=Ti[p(t)Oj], j = 0,l,...,L, (18) 

(Oo = I is the identity operator) and the normalization in order to maximize the entropy, 
defined (in units of the Boltzmann constant) by 

S(p) = -Tr[p\np). (19) 



Eq. |17| is a generalization of the more familiar density operator. For e.g. in open system, 
where we have Grand Canonical Ensemble there are two Lagrange multipliers, (3 = and 
/i are present, and we write the density operator as [^(J 



p(t) = exp (j3Q{T, V, p)-(3H + 0nN) , 



(20) 



As is well-known the dynamics are governed by the time evolution of the statistical operator. 
The time evolution of the statistical operator is given by 

i%% =[H(t),p(t)\. (21) 



The essence of the MEP formalism in conjunction with the group theory method is to 
find the relevant operators entering Eq. [17]) so as to guarantee not only that S is maximum, 
but also is a constant of motion. Introducing the natural logarithm of Eq. [17] into Eq. |21]) it 
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can be easily verified that the relevant operators are those that close a semi-Lie algebra 
under commutation with the Hamiltonian H, i.e. 

[H{t\6 j ]=inY l g ij {t)6 i . (22) 

i=0 

Thus the relevant operators may be defined as those satisfying the above equation. Equa- 
tion ^) defines an L x L matrix G and constitutes the central requirement to be fulfilled 



by the operators entering in the density matrix. The Liouville Eq. |21] can be replaced by 



a set of coupled equations for the mean values of the relevant operators or the Lagrange 



multipliers as follows [21 



d<0 j >t = -jr 9iS <d i >, j = o, i, (23) 

dX L 

-^ = E A ^> j = 0, 1,...,L. (24) 



In the MEP formalism, the mean value of the operators and the Lagrange multipliers belongs 
to dual spaces which are related by [|H 



IV. THE RELEVANT OPERATORS AND EVOLUTION EQUATIONS FOR THE 
HUBBARD HAMILTONIAN WITH A MAGNETIC FIELD TERM 

For notational convenience we now drop the subscript i, in all formulae from now on. 
The set of relevant operators for the HH with a magnetic field term is more than twice the 
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number of relevant operators without it. It is thus informative and useful to give the set of 
the relevant operators for the HH without the magnetic field. To this end we first consider 
the Hamiltonian form given in Eq. |12|. A little work shows that number operator n does not 
commute with the Hamiltonian, after some calculation we obtain 



-irp, 



(26) 



where p is given by 

p d ^ i(A*c\ + AJc| - A T c T - A&). (27) 

p is the mean field electron's current. Thus so far we have introduced three operators besides 
the Hamiltonian, namely h, x and p belonging to the relevant operator set. To determine 
the whole set we must proceed by finding the commutation relations of all the operators 
with the Hamiltonian until we get the complete set. The commutation relation of x with 
the Hamiltonian yields 

H,x] = -iap-iUL, (28) 



where /_ is the mean field pair's current and can be written as 

L = i([A*c\ - A T c T ]^ + n T [AJcj - A^q]). (29) 

We note that since the c's are fermion operators they anticommute, hence as a consequence 
of this 6j commutes with h^, viz, explicitly, c^n-f = Cjc|c| = — cjqcf = cjcfq = n^q. 

The commutation relation of p with the Hamiltonian introduces yet another two op- 
erators l + and Cj\. l + represents the mean field pair's interaction. 
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H,p = iax + iUl + — i4TUi, 



(30) 



where l + and 6)i are defined as 



l + d ^ f ([A*c] + A T c T ]% + n T [AJ£l + A 4 6J), 



•at 



UJ-y 



dcf |^tl 



c T ,c{ 



+ 



|AJ 2 



+ A^A^c} + A T AJc T c{ 



(31) 



The commutator of /_ with the Hamiltonian yields the final relevant operator of the 
present set, namely a) 2 > 



H, 1- = i(a + U)l + — iru)2, 



(32) 



d>2 is given by the following expression 



„ def „, 

u; 2 = 2( 



|Al 



c T ,cj 



"I + 



|A,I 2 



-n-r 



A^c^ + A T AJc T c{, 



(33) 



The three remaining commutation relations required to close the algebra can be expressed 
entirely in terms of operators already defined. These read 





= -i(a + U)L 




= i2\A\ 2 Tp, 


H,LJ 2 


= z8|A|VL. 



(34) 



Thus we have a set of seven relevant operators, namely h, x, p, /_, l + , u>i and u> 2 



which close the algebra as is clear from Eqs. go[ |28], |3(| [32| and |3J. Using Eqs. £2] and £3[ 
the evolution equations for the present set of relevant operators immediately follow and are 
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d <n> t 
dt 

d <x> t 
~dt 

d <p> t 
dt 

d < 1+ > t 
dt 

d < L > t 
dt 

d <LJi> t 
dt 

d <uj 2 > t 
dt 



]i <P>t > 
a 



U 



J- <P>t +JT <1->U 

a XJ ? 4r 

-- <x> t <l+> t +-z- <Wi> t , 



{a + U) 



n 

{a+U) 



<l->t, 



< l ->t +T <U 2 >t, 

n a 



2|A|V 
8|A| 2 r 



<L>t, 



(35) 
(36) 
(37) 

(38) 
(39) 
(40) 
(41) 
(42) 



The magnetic field term modifies the HH by a simple looking term, viz N, as is 
immediately apparent from the Hamiltonian form in Eq. We observe that N differs by a 
negative sign between the number operators of spin-up and spin-down states from h. This 
observation leads us to expect that like n, N when commuted with the Hamiltonian will 
lead to a set of relevant operators parallel to the ones obtained in case of h. 



H,N 



-irP, 



(43) 



where P is given by 



P ^ *(A*c{ - AJcI - A-fCf + A ; cJ 



H,X = -iuP -iUL- -ihp + iArVLx 



(44) 
(45) 
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where X, L_, Qi are given by 



L_ ^ i([A$ - A T c T ]^ - n T [AJc{ - A^cJ), 

fii d = z(A*AJc{c{ + AtA^cj. + A T AJc T c{ + AJA^cJ. 



iaX + iUL+ + ihx — i4rfl2, 



where L + , Cl 2 read 



L+ * ([A*c{ + Atotjni - n t [AJc{ + A^J), 



fi 2 = 



def , |A|| 2 



+ 



|AJ 2 



in 



+ A*AJc{cI + A i A T qc T ). 



The commutators of L+ and L_ with the Hamiltonian are 



H, L + =-i(a + U)L^-ihL-i2'rCl 1 , 



and 



if, L. 



i(a + £/)L+ + ihl+ + i2rtt 2 - 



The commutators of n, x, p, with the Hamiltonian in presence of magnetic 



H,h 
if, x 
H,p 
Hj + 



= -irp, 

= —iap — iUl- — ihP, 

= iax + iUl + — iArui + ihX, 

= — i(a + U)l- — ihL_, 
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H, L 


= i(a + U)l+ — iru)2 + ihL + , 






= i2\A\ 2 Tp- i2h(l 3 , 




H,u 2 


= i8|A| 2 rL - iAh(l 3 . 


(51) 



If we set h to zero we recover the equations obtained before, which provides a check on our 
calculations. Q 3 and f2 4 are defined as 



Cl 3 d ^ Fi(AJA t c{c T + A i A|c 1 c" 
64 = AJA^-A^cf. 

In addition we define two more operators Q5 and Qq 

fi 6 & *(A*AJc|c| - A^c^ 



(52) 



(53) 



It is clear from the definitions of tii, Q 3 and f2 6 given respectively in K6], |52| and f3 that 



fi 6 = ^ + ft 3 . 



(54) 



Eq. |54] provides a check on our calculation since it implies that once we have independently 
calculated the time-evolution equations for Cl 6 , Cli, and ^3 they must obey the relation 



d <Q 6 > t d <Qi> t d <Q 3 > t 



dt 



dt 



dt 



(55) 



Similarly it follows from definitions of Cl 2 [see f48j, S7 5 [see |53[ and the definitions of n and 
N that 



Q, = Q 2 + \A\ 2 I - \A\ 2 n - \A\ 2 N 



(56) 



which implies that 
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d <VL 5 > t d <VL 2 > t 



dt 



dt 



l ~ l9 d<n> t t ~ l9 d<N> t 
A - A 



dt 



dt 



(57) 



The commutators of f2i through f2 6 with the Hamiltonian are obtained after some 
calculation and may be displayed as 



h , 



H,Qt 



i(2a + U)Q 5 - i 2 r[| A| 2 x - |A| 2 X] + i 2 /i 4 , 

-2(2a + f/)fi 6 -2zr|A| 2 P, 

zr|A| 2 .x - zr|A| 2 X - z2/ift 4 , 

zr|A| 2 p-zr|A| 2 P-i2/ifi 3 , 

-z(2a + E/)fi 6 - ir|A| 2 P - zr|A| 2 p, 

z(2a + t/")fi 5 - ir| A| 2 £ + zr|A| 2 X. 



(58) 



It follows from the above discussion that we have a set of eighteen relevant operators 
in the presence of the external magnetic field, namely n, x, p, /_, l + , Cj\ , a) 2; N, X, P, 
L_, L + , Cli , ^2, ^3, ^4, ^5 , and which close the algebra as is clear from Eqs. BB, pj 



32], |34|, |43|, 03, 0, |50], |49] and |58|. However not all of the operators are independent as is clear 



from the relations given in Eqs. [53] and [56]. Using Eqs. |22] and EJ| the evolution equations 
for the present set of relevant operators can be written as 

d <n> t 



dt 

d <x> t 
dt 

d <p> t 
dt 

d < 1+ > t 
dt 



ol U 7 h " 

- <p> t + T <l-> t +T <P>t, 

nan 



-- <x> t 
n 



U 

h 



<l+>t 



At 



<UJ 1 > t 



h 



<x> t , 



[a + U) f h ± 
r <l- >t +t <L-> t , 
n n 



(59) 
(60) 
(61) 

(62) 
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d<l-> t (a + U) t r h ? . . 

3T- g = ~ t <l->t+T <& 2 >t-T <L+>t, 63 

at n n n 

d<u 1 > t 2\A\ 2 t h A 

^r^ = -^ <p> ' +2 A <n3> " < 64) 

d<X > t a * U $ h ta /n 

3T— " = T <P>t +t +t <p>< -4- <fii> 4 , 67 

at n n n n 

d<P> t a. U t h t A 

37— " = "T < X> t ~T < L +>t ~T <*>t +4- <V 2 > t , (68) 

at n n n n 

d<L + > t (a + U) t h f t ~ 

17^= z <L-> t +t <l-> t +2 ¥ <Vi>t, (69) 

at n n n 

-^ = -^<L + >«--</ + > 1 -2-< Si2 > t , (70) 
= ~ — j-^- <n 5 > t -2- <Q i>t +2-[|A| 2 <x> t -|A| 2 <X> t ], (71) 

= --|A| 2 <x> t +-|A| 2 <X> t +2- <fi 4 >t, (73) 

at n n a 

= --|A| 2 <p>i + -|A| 2 <P> t +2- <fi 3 > t , (74) 

d<n 5 >i (2a + £7) ^ r.r.n A r IAl2 - 

< Or >* (2a + U) ^ 7", r l9 r, . „ , , 

(77) 



To check the identity given in E3I we add Eqs. I7TI and 1731 and see if the sum of these two 
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equations agrees with Eq. [76]. We can immediately see that indeed this is the case. Likewise 
using Eqs. |5£], ^6], |72] and [75] we can see that the relation given in Eq. ^ holds. Thus we 
have an independent check of our stated relations. 

The following remark is in order in context of future outlook. The MEP formalism is 
limited to the mean-field approach. However the group theory based approach of identifying 
the set of operators which close the partial Lie algebra under commutation with the Hamil- 
tonian, is quite general. It is thus tempting to go beyond the mean-field formalism and use 
the set of relevant operators and their evolution equations to develop a technique which can 
take into account the quantum fluctuations. 

V. CONCLUSIONS 

We have given the set of relevant operators and the corresponding temporal evolution 
equations for the Hubbard Hamiltonian in the mean field approximation in the context of the 
maximal entropy formalism. The mean field approximation has been applied to the hopping 
term in the Hubbard term. As intuitively expected the inclusion of the external magnetic 
field leads to much larger set of the relevant operators than present in the case where the 
magnetic field is absent. 
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